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Abstract 

A theoretical model of stabilization of a microtubule assembly due to microtubule-associated- 
proteins(MAP) is presented. MAPs are assumed to bind to the microtubule filaments, thus prevent- 
ing their disintegration following hydrolysis and enhancing further polymerization. Using mean- 
field rate equations and explicit numerical simulations, we show that the density of MAP (number 
of MAP per tubulin in the microtubule) has to exceed a critical value pc to stabilize the structure 
against de-polymerization. At lower densities p < pc, the microtubule population consists mostly 
of short polymers with exponentially decaying length distribution, whereas at p > pc the average 
length increases linearly with time and the microtubules ultimately extend to the cell boundary. 
Using experimentally measured values of various parameters, the critical ratio of MAP to tubulin 
required for unlimited growth is seen to be of the order of 1:100 or even smaller. 
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I. INTRODUCTION 



Microtubules are hollow, cylindrical polymers which perform several important functions 
inside the cell. The persistence length of a typical microtubule is of the order of mm [ij], so 
they behave as rigid rods inside cells (linear dimension ~ 1 — lO/xm). The basic monomer 
unit of a microtubule is a heterodimer of a and /^-tubulin. Within a microtubule, tubulin 
heterodimers are arranged head-to- tail to form linear protofilaments, usually 13 per tube. 
The length of a tubulin dimer is nearly 8 nm, and the total diameter of a microtubule is 
about 25 nm. Heterodimers are oriented with their /3-tubulin monomer pointing toward the 
plus end, or fast growing end, of the microtubule. The minus ends are at the microtubule- 
organizing centers such as centrosomes in animal cells. 

Both a and /5-tubulin binds GTP. Following polymerization, the /3-tubulin in the dimer 
hydrolyses its bound GTP to GDP |2| . The hydrolyzed GDP-tubulin (which we shall refer to 
as D-tubulin, while the un-hydrolyzed version shall be called T-tubulin) does not polymer- 
ize, and the process is irreversible inside the microtubule. When the advancing hydrolysis 
front reaches the growing end of a microtubule, it starts de-polymerizing from the end and 
the protofilaments start peeling off, releasing the GDP-tubulin into solution. Outside a 
microtubule, reverse hydrolysis takes place and polymerization events start all over again. 
Thus the microtubule constantly switches between phases of growth and shrinkage, and this 
behavior (first observed in in vitro experiments) is called dynamical instability^]^. When 
a microtubule makes a transition from growing to shrinking state, it is said to undergo 
catastrophe and the reverse transition is called rescue. The relative rates of catastrophe and 
rescue, combined with the velocities of growth and shrinkage determine the character of a 
given population of microtubules j^O,!?!- 

The intrinsic microtubule dynamics observed in vitro is modified inside the cell by interac- 
tion with cellular factors that stabilize or destabilize microtubules, and which operate in both 
spatially and temporally specific ways to generate different types of microtubule assemblies 
during the cell cycle. Proteins that modulate microtubule dynamics are called microtubule- 
associated proteins (MAP). Several MAPs have been identified, including MAPI, MAP2 and 



tau in neurons and MAP4 in non- neuronal cells j8|, |9|, |lO|, |ll| . MAPs regulate the microtubule 
dynamics through one or more of the following mechanisms: suppressing catastrophe rate, 
increasing rescue rate, promoting polymerization and preventing de-polymerization 
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Recent experimental studies show that, inside cells, dynamic instability of microtubules is 
observed predominantly near the cell margin. Deep inside the cytoplasm, microtubules were 
observed to be in a state of almost persistent growth with very few catastrophe events. Direct 
experimental measurements showed a marked difference between the rates of catastrophe 
in the cell interior (z/c ~ 0.005s~^) and periphery (z/^ ~ 0.08s~^)[13]. These observations 
suggest that in the cell interior, microtubule growth is sustained by cellular factors which 
are presumably absent or rare near the cell boundary. Our principal aim in this paper is 
to develop a simple mathematical model to study how the microtubule-associated proteins 
affect the dynamics of microtubule assembly. 

It is convenient to outline the main results of this paper at this point. By analyzing the 
rate equations associated with microtubule dynamics under a mean-field treatment of MAP 
density, we show that the microtubule assembly can be in two different phases depending 
on the proportion of MAP to tubuhn (henceforth referred to as the MAP density) inside 
the microtubule. Above a critical density pc, individual microtubules exist in a phase of 
unlimited growth, while below this density, growth is limited and the average length is 
finite. We show that pc = Pu/iPg +Ph), where pn and pg are rates of hydrolysis and growth 
respectively. The rates of catastrophe and rescue arc related to p through the relations 
J^c — Ph{^ — p) and Ur — PsP where Ps is the rate of shrinkage. Numerical simulations show 
that the fluctuations in MAP distribution do not signflcantly affect the predictions based 
on mean-field analysis. We compare our results to experimental observations on the effects 
of MAP on microtubules, and find that there is qualitative agreement. 

This paper is arranged as follows. In the next section, we briefly discuss the salient 
features of our model and its limitations. In Sec. 3, we write down explicitly the rate equations 
that govern the dynamics of the microtubule assembly. We derive expressions for the critical 
MAP density, the time evolution of the average length of microtubules in the unlimited 
growth regime as well as the expressions for the rates of catastrophe and rescue in our 
model. In Sec.4, we present results of numerical simulations of our model which are in 
excellent agreement with the analytical predictions. In Sec. 5, we summarise our results and 
outline the main conclusions that emerges out of this study. 
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II. THE MODEL 



In this model, we assume that microtubules are thin rigid rods nucleating from the 
microtubule-organizing center, which is very often the centrosome(Fig.l). The three- 
dimensional structure of microtubules is ignored in our model. The centrosome contains 
nucleation sites for a microtubule, and nucleation takes place at vacant sites at a rate u. 
The microtubule grows by adding tubulin sub units at a rate p^, and tubulin units inside 
the microtubule are stochastically hydrolyzed at a rate pn- A microtubule will start de- 
polymerizing when the T-tubulin at its growing end is hydrolyzed to D-tubulin ^, 14 , Q| . 
The shrinkage takes place at a rate Ps-, which is the probability per unit time to lose a 
subunit. It is generally assumed that the de-polymerization process will continue along the 
length of the microtubule, until the shrinking end encounters a patch of T-tubulin in the 
interior. In our model, however, we assume that a microtubule in shortening state will con- 
tinue to shrink (irrespective of the state of tubulin in the interior) until the de-polymerizing 
end encounters a MAP, which are present along the microtubule. This simplification helps 
focus on the role of MAPs in the microtubule dynamics. We further assume that MAPs also 
assist in the polymerization process by acting as polymerization sites along a microtubule. 
For further simplification, we assume the presence of an excess of free tubulin in solution, so 
that the dynamics of the tubulin concentration can be largely ignored. We also assume that 
the nucleation sites in the centrosome are spaced sufficiently far apart so that the effective 
interaction between microtubules through local depletion and enhancement of concentration 
of free tubulin is insignificant. A schematic illustration of the dynamics of our model is given 
in Fig. 2, and a schematic illustration of the various rates of the model is shown in Fig.3. 

It is helpful to note the typical experimental values of the various parameters in our model. 
We refer to recent experimental work [13] for most of the parameter values. Nucleation of 
microtubules in the centrosome was seen to occur at a rate of ^ 5-6/min in vivo jisl ]. 
The velocity of growth of growing microtubules observed is Vg ~ 20fim/min. The rate 
of addition of tubulin dimers (pg) is determined from Vg as follows. A tubulin dimer has 
length ^ 8nm. Since there are 13 protofilaments of T-dimers in a single microtubule, each 
tubulin contributes a length of 5 ~ 0.6nm to the microtubule. Thus, the rate of addition 
of tubulin dimers would be Pg = Vg/6 ~ 0.55 x lO^s^^. Using similar reasoning, the rate 
of loss of tubulin can be calculated from the measured velocity of shrinkage in shortening 
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microtubules. This turns out be ps = Vg/S ~ 0.82 x 10^s~^. 

The rate of hydrolysis is a difficult parameter to measure, and is not as well known as the 
other rates. However, a wide range of values, 0.25-25/min has been reported 
The typical values of the density p of MAP used in injection experiments are estimated to 
be of the order lO'^ WmJ. 



III. MEAN-FIELD RATE EQUATIONS 

In this section, we analyze the mean-field rate equations for the time evolution of the 
probability distribution of the lengths of microtubules. Since our basic aim is to study 
the character of an assembly of microtubules, we resort to a statistical description. The 
most appropriate quantity to study in this perspective is the length distribution p{l,t), 
which is the fraction of nucleation sites which has a microtubule of length / at time t. The 
presence of MAP is taken into account through a mean-field density variable p, whose local 
fluctuations in space is ignored (The effect of these fluctuations are explicitly studied in 
numerical simulations presented in Sec. 4). The MAPs change the microtubule dynamics 
by modifying the rates of shrinkage and hydrolysis, as seen below. It is also convenient 
to classify the entire set of microtubules into two: growing and shrinking microtubules (as 
was done in [151]). Consequently, the complete length distribution p{l,t) can also be split as 
follows. 



p{i,t)=p4i,t)+p4i,t) (1) 

where p+{l,t) is the fraction of microtubule in growing state at any time t, and p_{l,t) is 
the fraction in shrinking state at time t. The rate of change of p+{l,t) is described by the 
following equation (for / > 2): 



^P±^=p_^{l-l,t)p,+p4l + l,t)ppg-p4l,t)[pg+pn{l-p)] ; l>2 (2) 

The significance of each term easily follows from the model description in Sec. 2. The first 
two terms represent gain events, where a microtubule of length / in growing state is created 
at time t. In particular, the second term describes the transition between a microtubule in 
a shrinking state to a growing state, and so 



Ur = PPs (3) 

is the rate of rescue in our model. We note that in our model z/.,. ^ as p — > 0, which is 
in agreement with the assumptions of our model. (However, in experiments, microtubules 
have been observed to undergo rescue under in vitro conditions in the absence of MAP. Our 
assumption here is that the MAP significantly increases the rescue frequency, which would 
otherwise be negligible). The last two terms represent the loss events where a microtubule 
of length / either adds a tubulin to increase its length, or undergoes hydrolysis to transform 
into shrinking state. The latter is an event of catastrophe, and so the corresponding rate in 
our model is 

u, = ph{1-p). (4) 

Both the rates of catastrophe and rescue are parameters which are measured in experiments. 
We also note that for microtubules of unit length, the gain event in p^{l,t) is nucleation, 
rather than polymerization, so that its equation is 

^^^^ = ^Rit) + P-(2, t)pps - t) [p, + pj,{l - p)] (5) 

where R{t) = 1 — Y^'^i[p+{l,t) is the fraction of vacant sites in the lattice at time 

t. The equation for p^{l,t) is, 

^P^^=p^{l,t)pH{i-p)+p-{l + ht)Psii-p)-P-{l,t)Ps ; l>i (6) 

The equations Eq.(j21) and Eq.® are supplemented by the boundary condition for p{l,t), 
which is, 

limi^ooP{l,t) = (7) 

We now proceed to look for the solution for these coupled set of rate equations. As a first 
step, we look for the steady state solution by putting the time derivatives to zero in Eq.Q 
and Eq.dni). 
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A. The steady state: 



In the steady state, it is reasonable to assume that both the distributions p+{l) and P-{1) 
(the time dependence has been dropped since we are looking at steady state solutions now) 
have the same dependence on /. Further, it is easily seen from the equations that the steady 
state solution has an exponential form, which is easily determined by substituting the trial 
functions 



p+{l) = Ae-'^ ■ p_{l)=^p^{l) (8) 

in the steady state equations. In the trial solution, A and 7 are two unknown constants which 
need to be determined, and A is the normalization factor. Upon substitution in Eq.(jHI) we 
find the following equations for the unknown constants A and 7. 

Pg[e^ - 1] + ipspe-^ = Ph{1 -p)= Psl[l - (1 - pje"^)] 
from which we determine the unknown constants A and 7 as follows. 

A = %[(l + ^)(l-p)] ; 7 = ^^±^(l-p) (9) 

Pg Ps 

Clearly, from the boundary condition in Eq©, the solution is meaningful only if A > 0, and 
this implies that p < pc, where 

Pc = —— (10) 

Pg+PH 

It turns out that pc is a critical value of MAP density, below which growth is bounded and 
the length distribution of microtubule is stationary. The typical values of pc can be easily 
estimated using the known values of pn and Pg. For example, keeping pg ~ 0.55 x lO^s^^ 
and changing pn from ~ 3/min to 300/min (which are close to in vivo estimates) gives pc 
in the interval ~ 10~^ — 10~^. This suggests that a very small MAP to tubulin ratio might 
be enough to stabilize microtubules against depolymerization. 

The steady state solution discussed above is not valid when p > pc, and so the solution 
for the rate equations has to be found explicitly for this case. 
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B. The unbounded growth regime 



When p > pc growth becomes unbounded. Although the exact solution of the evolution 
equations in this regime is hard to find, it is helpful to note that under usual experimental 
conditions the growth and shrinkage rates far exceed the hydrolysis rate, i.e.,pH is small in 
comparison with pg and Ps- If we further assume that we work in the regime where the 
MAP density is high, i.e., p ^ 1, then all the terms in the rate equations which contain 
Ph can be neglected in comparison with the other terms. Further analytic progress is now 
possible by converting the finite difference equations into differential equations (in the limit 
/ ^ 1). This is done by formally expanding t) in its derivatives (see the appendix for a 
discussion on this point) as the series, p±{l±l,t) = pj-{l,t)±^^ + ^^^^. After substitution, 
we arrive at the following differential equation for p_{l,t). 

dp_{l,t) Ps d'^p-{l,t) .dp-{l,t) 

-^r- = Y^' - - p)^^-p-(^^')p^p (11) 

Let us now define the re-scaled time r = ppst and the new variables D = v = 

in terms of which Eq. (fTT|l is re-cast as 

dp-{l,T) d^p4l,T) dp^jl.r) 
^^ = ^^^ + "^^ 

This equation is easily solved by substituting p_(/, r) = e^'^g{l, r), upon which we find that 

g{l, r) satisfies the equation |^ = + v^, which is the well-known equation for diffusion 

of a particle in one dimension, with a drift term. The solution of this equation (apart from 

the trivial solution g{l, r) = 0) has the form(up to a proportionality factor), 

1 (I + vt)^ 

g{l,r) = ^=e-^i^ (13) 

However, this solution is not physically reasonable since the mean length decreases with 
time as (/) = —vt. It follows that in this limit, the solution for Eq. ()ll|) is the trivial solution 
i.e, p^{l,t) ~ 0. The equation for p^{l,t) in the same limit is, 

dp+{l,T) _ Pg d'^p+{l,T) dp+{l,T) 

dt 2 dP dl ^ ' 

which, again, is the diffusion equation with a drift term, whose solution is (up to a constant) 
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p+il.t) = 




We see that in the hmits pn ^ Ps and p ^ 1, the length distribution of growing mi- 
crotubules is a Gaussian at all times, and the mean value increases linearly with time as 
(/) ~ Pgt. Furthermore, in this limit, almost all the microtubules are in the growing state. 

Although it is possible at this stage to look for a perturbative solution of the evolu- 
tion equations Eq.()ll|) and Eq.()14|) with ^ as the pertrubation parameter, we have not 
attempted this here. Rather, in the next subsection, we determine the mean length as a 
function of time exactly, and shows that it agrees with the prediction of Eq. p5|l in the limit 
Ph < Ps- 



C. The average length of microtubules 

In order to study the time evolution of the average length in the unbounded growth 
regime, it is convenient to define the following quantities: 

oo 

g+(t) = 5:p+(/,t) (16) 

1=1 

which is the fraction of microtubule in growing state, irrespective of length, and 

oo 

Q^{t) = Y.p_{l,t) (17) 
1=1 

is the fraction of microtubule in shrinking state, irrespective of length. The time evolution 
of Q+{t) can be found by performing a sum of Eq. ^ over all / > 2 and adding Eq.(j3)).The 
result is 

= uRit)~pHil-p)Q+it)+pps[Q4t)-p4l,t)] (18) 
For Q-{t), we perform a sum over Eq.® over all / > 1, and find that 

=p^(l -p)Q+(t) +p,(l -p)[Q_(t) -p_(l,t)] -psQ-it) (19) 
where R{t) = 1 — Q+{t) — Q^{t). The average length of microtubule at time t is defined as 
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{i){t) = Y.i[p4i,t)+p-ii,t)] (20) 
1=1 

The rate of change of (/) with time may be expressed in terms of Q+{t) and Q^{t), since 
each growing microtubule grows by an amount Pg and each shrinking microtubule shrinks 
by an amount Ps per unit time. It follows that 

^ = Q+it)pg-Q_it)p, (21) 

To calculate Q+(t) and Q~{t) from Eq. (|THjl and Eg. ([11^1 . let us first assume that t) <^ 
Q-{t). This assumption is valid in the unbounded growth regime, where the mean length 
increases with time, and hence the peak of the distribution is expected to shift toward larger 
and larger times (note the shifting Gaussian we found in the last section). The result is a 
closed set of equations for Q+{t) and Q-(t), which has the form 

^^j^ = u{i - g+ - g_) - «g+ + (3q^ (22) 

and 

= aQ. - PQ^ (23) 

where a = Ph(1 — p) and /3 = pps- We take the initial condition to be Q+{0) = v and 
g_(0) = 0. The general solution for the above set of equations is 

and 

^ _ ^ , az/(l-(a + ^)) a{\-v){v-p) ^ 

^ « + (a + /?)(« + /?-z/) ^ (« + /3-z/) ^ ^ 

We note that the exponential terms vanish at sufficiently late times, and so the asymptotic 
forms of g+(t) and Q-{t) are. 



Q-it) 



PsP 



Psp + Ph{1 - p) 
^ Ph{1- p) 
PsP + Ph{1- p) 



t ^ 00 
t ^ 00 



(26) 
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Thus, even in the unbounded growth phase, a finite fraction of microtubule will be in the 
shrinking state. After substituting Eq. ^Bj) in Eq.(j2H), we find that 



dt pPs+PhO-- p) 
For unbounded growth, we need ^ > 0, for which the condition p > pc, where pc is 
defined in Eq. fllOj) . This calculation confirms that pc is the critical MAP density above 
which microtubules exist in a phase of sustained growth. In the limit pn Ps, we see that 
(Z) ~ pgt, which is one of the results we had derived in the last subsection. 



D. Comparison with experiments 



At this stage, it is worthwhile to compare our predictions (although based on a very siraple 
model) with experimental observations regarding the effect of MAPs on microtubules [isl. Ifoj] . 
We specifically refer to experiments by Dhamodharan and Wadsworth jl^ where MAPs were 
microinjected into BSC-1 cells and their effects on growth of microtubules, in particular, the 
changes in catstrophe and rescue frequencies, were studied. When 0.1 mg/ml of MAP-2 
was added to the cells, the proportion of MAP to tubulin in the polymerized state was 
estimated to be nearly 1:25, or p ^ 0.04. Increasing the MAP concentration 10-fold was 
observed to decrease the catastrophe frequency from 0.05 (±0.02)s^^ to 0.02 (±0.003)s^^, 
i.e., by a factor of nearly 0.4 after neglecting the experimental error. According to Eq. (j^, 
increasing p by a factor of 10, from 0.04 to 0.4, would decrease the catastrophe frequency 
by a factor of 0.625, which is in resonably good agreement with the experimental value. 
Furthermore, according to Eq. the catastrophe rate in the absence of MAP is simply 
the hydrolysis rate pn, which may be calculated from the same equation using experimental 
data for Uc and p. The data for p = 0.4 (which has greater accuracy) gives pn ~ 0.033s~^, 
which is, again, in reasonably good agreement with the experimental value of 0.043 (±0.01) 
s~^. The effects of MAP on the rescue frequency is more subtle, since rescue would occur 
in cells even in the absence of MAP, which is not taken into account in our model. 
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IV. NUMERICAL SIMULATIONS 



In this section, we discuss the results of our numerical simulations of the microtubule 
assembly. Since our model neglects any effective interaction between growth of neighboring 
microtubules due to local depletion of hgand concentration, the geometry of the system, and 
the spacing between individual microtubules is unimportant in the simulation. For simplic- 
ity, the cell interior is imagined as a rectangular box, where the base forms the substrate 
for nucleation of microtubules |15|]. The microtubules grow in the direction perpendicular 
to the base until they hit the boundary. 

The next step in the simulation procedure is to divide the box into a grid array of 
dimensions L x L x H, where L is the length of the side of the substrate and H is the 
height. As we had mentioned before, the lattice spacing in the substrate is unimportant as 
long as it is sufficiently large that direct interaction between neighboring microtubules could 
be ignored (henceforth, we shall always assume that this condition is satisfied). The lattice 
spacing in the direction of growth (z-axis) is the length of one tubulin dimer, which, after 
taking into account the number of protofilaments, is equivalent to 6z ~ 0.6nm. For the time 
evolution, we choose Monte Carlo step as equivalent to 6t ~ 10~'^s. With this timescale, 
the probability of a microtubule growing by one unit per MC step is nearly 0.55, and the 
probability that a tubulin unit is lost in one MC step is nearly 0.82, from the rates discussed 
in Sec. 2. In other words, the time scale has been chosen such that growth and shrinkage 
of microtubules are stochastic events, in accordance with the rate equations. For the lattice 
size, we chose L = 100 and H = 500. 

The tubulin density in polymerized state at each point in space is represented using a 
density index 0(r, z, t), where r is the substrate co-ordinate, z is the position in the direction 
of growth and t is time in units of MC steps. At t = 0, we start with a substrate free of 
microtubule, i.e., 0(r, z, 0) = 0. We do not explicitly keep track of the density of free tubulin 
in the solution. The density index (j){r,z,t) takes three values: (f){r,z,t) = +1 represents 
T-tubulin, (j){r,z,t) = —1 represents D-tubulin and (f){r,z,t) = implies that the site is 
empty. A string of (f){r,z,t) 7^ at any r, for 1 < z < / is a microtubule of length /. For 
z = 1, 0(r, z, t) changes from to +1 with probability u per MC step at all r and t. This is 
the nucleation event. For any z > 1, this transition takes place with probability pg, which 
is a polymerization event where a new T-tubulin dimer is added to the microtubule. The 
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location of MAP in the lattice is specified through another density index x(r, t), which takes 
the value 1 or depending on whether there is a MAP molecule at the lattice site (r,z) at 
time t. At all times, we distribute the MAPs at random throughout the lattice with mean 
density p. We do not take into account the diffusion of MAPs in solution explicitly. Rather, 
when microtubules grow from the substrate, they assimilate the MAPs along their length 
and release them randomly such that the mean density remains p at all times. We assume 
here that the binding and release of MAP take place over time scales much less than one 
IMC step in the simulation (which is equivalent to 10"'^ seconds), so that, at each time step, 
the microtubule assembly has a different MAP configuration. 

In all our simulations, the rates of growth and shrinkage were fixed as Pg = 0.55 and 
Ps = 0.82, which are the experimental values quoted in Sec. 2. For reasons we shall shortly 
discuss, the hydrolysis rate was fixed at a value higher than the typical estimates, i.e, we 
chose ph = 0.5 (which would translate into a hydrolysis rate of ~ 3 x 10^/minute). We also 
chose a high nucleation rate (z^ = 1), for two reasons: to increase the occupancy of the lattice 
so as to improve the statistics, and to reduce the duration of the early-time regime of growth 
(where the exponential factors in Eq.. ()24|l and Eq.. ()25j) would ne non-negligible). The chosen 
values of pn and Pg gives a convenient value of the critical MAP density, p^ — 0.47 (which 
was the primary reason for choosing a high pn)- 

We measured the length distribution and average length of microtubule assembly for three 
values of p, 0.1, 0.3 and 0.6. We chose the first two values to be below pc, where we expect 
bounded growth (steady state distribution) and the last value is above pc, where growth is 
expected to be unbounded. The results for probability distribution of lengths are given in 
Figs. 4-6, and the result for time evolution of the average length in the unbounded regime 
is displayed in Fig. 7. In Fig. 4 and 5, we observe that, as predicted by the rate equation 
analysis, the solution is exponential, but the measured value of A differs from the mean- 
field prediction. This is presumably due to the spatial fluctuations in the MAP density, 
which is ignored in the mean-field analysis. In Fig. 6, we find that the length distribution of 
microtubules in the unbounded regime closely resembles a Gaussian function with a mean 
value shifting with time. This is in accordance with the results of rate equation analysis in 
Sec. 3 (although in the simulations pn and Ps are of the same order of magnitude). In Fig. 7, 
we confirm that the mean length in the unbounded regime follows a linear dependence on 
time, again in agreement with the analysis in Sec. 3. 
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V. SUMMARY AND CONCLUSIONS 



In this work, we have presented a simple theoretical model of stabilization of microtubules 
due to microtubule-associated proteins. We have shown that when the density of MAP is 
reduced below a critical value, the length distribution of the microtubule assembly under- 
goes a significant change. Whereas in one phase (at high MAP density) the microtubules 
extend throughout the cell, in the other phase they have an exponentially decaying length 
distribution with a finite average length. It is interesting to note that this transition is sim- 
ilar to the transition at mitosis in real cells, where phosphorylation of MAPs by agents such 
as the metaphase promoting factor (MPF) leads to a weakening of their binding to the 
microtubules, and a consequent change in length distribution. We have also shown that the 
critical ratio of MAP to polymerized tubulin required for a transition to the 'cell-spanning' 
phase can be as low as 1:1000 under in vivo conditions. Thus, our conclusion is that even a 
very small concentration of MAP can sustain unlimited growth of microtubules, which is in 
qualitative agreement with experiments jisl . 

There are still many open questions regarding the effect of MAPs on microtubule dynam- 
ics in cells. In fact, MAPs which inhibit polymerization have also been identified recently 



It would be interesting to include such MAPs also in our model and look for the conse- 
quences of competition between the effects of two kinds of MAP. Moreover, our prediction 
of the existence of a critical MAP density which differentiates between two distinct growth 
regimes for microtubules could be verified in future experiments. Lastly, our model could 
be still improved by including effects of local fluctuations in the concentrations of MAP and 
tubulin. In particular, it would be interesting to ask in the context of such an improved 
model whether there exists subpopulations of microtubules in one growth regime, while the 
over-all growth behavior of the whole population is in the other regime. These are some of 
the questions we would like to address in the future. 
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APPENDIX 



In this appendix, we outline the scheme by which the difference equations Eq..((21) and 
Eq..® are converted to partial differential equations Eq.. (fTT|) and Eq.. ()14|) . We first define 
the first derivative in a symmetric way by taking the average over right and left derivatives. 

dp± ^ [p±il + 1) - P±{1)] + [p±il) -P±{1- 1)] 
dl 2 

We note that the second and third terms cancel each other, so we arrive at 

dp± _ P±{1 + 1)-P±{1-1) 
dl 2 

For the second derivative, we use the familiar Euler discretization scheme, i.e, 

^^p±(/ + l)+p±(/-l)-2p±(/) 
After solving these two equations together, we arrive at the expression 

M/±i,t)=MM)±^ + 2^ 
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FIG. 1: A schematic illustration of microtubules inside the cell. The filled circle is the microtubule- 
organizing center, which is very often the centrosome. 
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FIG. 2: A schematic illustration of the dynamics of our model. The dark squares are the T-tubulin 
units and the white squares are tlie hydrolysed D-tubuUn units. The large filled squares represent 
the microtubule-associated proteins. In (1), hydrolysis of the microtubule has started from the 
top end, which spreads to the bottom (2), and untimately leads to the D-tubulin units falling off 
the polymer. But the de-polymerization is arrested in (4) by the MAP, which promotes further 
polymerization. 
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FIG. 3: A schematic illustration of the various rates in the microtubule dynamics. Mirotubules 
make the transition between growing and shrinking states with rates which depend on the MAP 
density. 
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FIG. 4: Probability distribution of microtubule length when p = 0.1. We observe bounded growth 
with exponential distribution of lengths as predicted by the rate equations. The unit of time in 
this plot is one second. /(/) oc e^^' is a fit function where A ~ 0.79, which differs somewhat from 
the mean-field prediction of Eq..Q, which gives A ~ 0.54. 
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FIG. 5: The figure shows the probabihty distribution of microtubule lengths when p = 0.3. In this 
case, growth is again bounded, but the length scales reached by the microtubule is longer than 
before. The unit of time in this plot is one second. /(/) oc e~^' is a fit function, with A « 0.49. 
The mean-field prediction for A given by Eq .@ is nearly 0.28. 
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FIG. 6: The figure shows the probabiUty distribution of microtubule lengths when p = 0.6 at three 
different times. The growth is unbounded since p > pc- The peak of the distribution shifts to the 
right with increasing time, which shows that the average length increases with time. The unit of 
time in this figure is one second and length is measured in units oi 5 = 0.6nm. 
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FIG. 7: The plot shows the average length of the microtubule assembly increasing linearly with 
time for /9 = 0.6, which is above pc- The straight line is a linear fit. Time is measured in seconds 
and length is measured in units oi S = 0.6nm. 
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